Shear viscosity and Bose statistics: 
Capillary flow above A point 

Shun-ichiro Kobe 

Physics Division, Faculty of Education, Kochi University 
Akebono-cho, 2-5-1, Kochi, 780, Japan 
(Dated: February 1, 2008) 

Abstract 

The gradual fall of the shear viscosity at T\ < T < 2.8K, observed in a liquid helium 4 flowing 
through a capillary, is examined. The disappearance of the shear viscosity in a capillary flow is a 
manifestation of superfluidity in dissipative phenomena, the onset mechanism of which is a subtle 
problem compared to that of superfluidity in non-dissipative phenomena. Applying the linear- 
response theory to the reciprocal of the shear viscosity coefficient 1/ry, we relate these two types of 
superfluidity using the Kramers-Kronig relation. As a result, we obtain a formula describing the 
influence of Bose statistics on the kinematic shear viscosity in terms of the susceptibility, without 
getting involved in the details of the dissipation mechanism of a liquid. Compared to an ordinary 
liquid, a liquid helium 4 above T\ has a 10 -3 times smaller shear viscosity coefficient. Hence, 
although in the normal phase, it is already an anomalous liquid under the strong influence of 
Bose statistics. The coherent many-body wave function grows to an intermediate size between 
a macroscopic and a microscopic one, not as a thermal fluctuation but as a thermal equilibrium 
state. Beginning with bosons without the condensate, we make a perturbation calculation of its 
susceptibility with respect to the repulsive interaction. Using the above formula on the kinematic 
shear viscosity, we examine how, with decreasing temperature, the growth of the coherent wave 
function gradually suppresses the shear viscosity, and finally leads to a frictionless flow at the 
A-point. 

PACS numbers: 67.40.-w, 67.20. +k, 67.40.Hf, 66.20.+d 
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I. INTRODUCTION 



Superfluidity of a liquid helium 4 was first discovered in a flow through a narrow channel 
and a capillary [1]. In an ordinary liquid, the shear viscosity causes the shear stress F xy 
between two adjacent layers moving at different velocities 

F =n— (1) 

where 77 is the shear viscosity coefficient. In a narrow pipe with an inside radius a and a 
length L, the field of velocity driven by the pressure difference AP has a form such as 

. , a 2 — r 2 AP 

<-< r > = — — • < 2) 

where r is a radius in the cylindrical coordinates (Poiseuille flow). In the superfluid phase 
of a liquid helium 4, even when the pressure difference vanishes (AP = 0), one observes a 
non- vanishing flow (v z (r) 7^ 0), hence 77 = in Eq.(2). 

The shear viscosity of a liquid helium 4 has been subjected to considerable experimental 
and theoretical studies. At < T < T\, an oscillating or rotating probe immersed in a 
liquid helium 4 experiences the shear stress in a complicated manner (viscosity paradox). At 
temperatures well below T\, various excitations of a liquid are strictly suppressed except for 
phonons and rotons. Hence, it is natural to assume that phonons and rotons are responsible 
for the shear viscosity. By regarding these elementary excitations as a weakly interacting 
dilute Bose gas, the formula in the kinetic theory of gases well describes the shear viscosity 
of a liquid at T T x , which is a success of the quasi-particle picture in the many-body 
theory fl. 

For the shear viscosity at T > T\, however, the dilute-gas picture must not be applied, 
because the macroscopic condensate has not yet developed. Instead, one must tackle the 
problem of how Bose statistics affects excitations and dissipations of a liquid. Figure 1(a) 
shows 77 of a liquid helium 4 measured in a flow through a capillary at the vicinity of T\ 

nnn 

[3J[4][5J. One notes that 77 does not abruptly drops to zero at T\, but it gradually decreases 
with decreasing temperature from Tx + 0.7 K, in which the macroscopic condensate has not 
yet developed, and finally drops at the A-point. 

For the shear viscosity of a liquid, Maxwell obtained a simple formula 77 = Gt (Maxwell's 
relation), where G is the modulus of rigidity and r is a relaxation time, using a physical 
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FIG. 1: (a) The temperature dependence of the shear viscosity rj(T) of a liquid helium 4 measured 
in a flow through a capillary (after the data in Ref.3 and 5). (b) The temperature dependence of 
the specific heat C(T) of a liquid helium 4 (after Atkins, Liquid helium4, (Cambridge, 1959)). (c) 
C(T) of an ideal Bose gas with the same density as a liquid helium 4. 

argument j(| (see Appendix. A). In the short time scale, a liquid shows similar behaviors to 
that of a solid at the molecular level. In a liquid, G is determined mainly by the motion of 
vacancy as well as in a solid. Since no apparent structural transformation is observed in a 
liquid helium 4 at the vicinity of T\, G may be a constant at the first approximation. The 
gradual decrease of rj = Gr above T\ suggests that the influence of Bose statistics gives rise 
to a considerable decrease of r. To understand superfluidity occurring in such a dissipative 
flow, we must deal with the microscopic mechanism of the dissipations determining r. At 
first sight, it seems to be a hopeless attempt, because the dissipation itself is a complicated 
phenomenon allowing no simple description [7;]. In a superfluid helium 4, many features 
associated to Bose statistics are masked by the strongly interacting character of the liquid. 
But, from a different viewpoint, new physics must lie in the phenomena in which Bose 
statistics is deeply entwined with the structure and dynamics of a liquid. Hence, we must 
find a starting point for disentangling these complicated problems. For the shear viscosity, 
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FIG. 2: Two types of flows, (a) a flow through a capillary, (b) a rotational flow in a container. 
White circles represent an initial distribution of particles. The long arrows move white circles on a 
solid line to black circles. In both (a) and (b), permutation symmetry holds over the whole liquid 
at T < T\. The particles represented by black circles are reproduced by slight displacements by 
short arrows from the initial particles close to that particle. 

we find such a starting point, in which one can consider one problem separately from the 
other ones. 

The viscosity has a formally similar expression to other processes of diffusion, such as 
diffusion of matters and thermal conduction. The latter are the diffusion of a scalar quantity 
such as density p and temperature T, and are always dissipative processes, whereas the 
viscosity is the diffusion of a vector quantity, momentum p. Since the vector field mv(r) 
has various spatial patterns, the viscosity has a qualitatively different feature that is not 
seen in diffusion of matters and thermal conduction. Figure 2 schematically illustrates two 
types of flows: (a) a flow in a capillary, and (b) a rotational flow in a container. In fluid 
dynamics, the viscous dissipation in an incompressible fluid is estimated by the dissipation 
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function 

$(r) = 2?7 (^ij - ^ekkdijj , (3) 

(eij = (dvi/dxj + dvj/dxi)/2 is the shear velocity). Particles in a capillary flow like Fig. 2 (a) 
experiences thermal dissipation not only at the boundary, but also in a flow. (For Poiseuille 
flow Eq.(2), the dissipation function Eq.(3) is not zero at any r except for r = 0.) On the 
other hand, in a rotating cylindrical container like Fig.2(b), a liquid makes the rigid-body 
rotation owing to its viscosity, the velocity of which is proportional to the radius and the 
rotational velocity Q such as v^r) = ft x r . Except at the boundary to the wall, there 
is no frictional force within a liquid, and the flow is therefore a non-dissipative one. (For 
v <i>{ r ) — Eq.(3) is zero at any r.) When we discuss the onset mechanism of superfluidity 
in such a non-dissipative flow, we can take an advantage of this feature to avoid dealing with 
the dissipation mechanism. 

In this paper, we will consider the frictionless capillary flow of a liquid helium 4 by relating 
it with the nonclassical rotational flow using the Kramers-Kronig relation. (The former is 
superfluidity as a non-equilibrium process, whereas the latter is superfluidity as a thermal 
equilibrium one.) This method originated in studies of electron superconductivity. Just after 
the advent of BCS model, an attempt was made to relate the electrical conductivity (more 
precisely, the micro-wave absorption spectrum) with the penetration depth in the Meissner 
effect js] {9]. The former is a quantity in the dissipative phenomenon, and the latter is that 
in the non-dissipative one. Furthermore, the Meissner effect has a common mechanism with 



the nonclassical rotational flow of a Bose liquid [10| . Hence, in the relationship between the 
electrical conduction and the Meissner effect in superconductivity, one finds a parallelism to 
that between the capillary flow and the rotational flow in a liquid helium 4. By using this 
analogy, one can consider the influence of Bose statistics on the shear viscosity of a liquid 
without getting involved in details of the dissipation mechanism of a liquid. Furthermore, 
one can explore the relationship between two typical manifestations of superfluidity . 

Before carrying out this program, one must have a physical image on the gradual fall of 
i] at T\ < T < 2.8K observed in a liquid helium 4 (Fig. 1(a)). It has been assumed to be 
a manifestation of the thermal fluctuation of Boson systems, and discussed in connection 
with the A-shape of the specific heat C( T) ( Fig. 1(b)). The principal peak in C(T) reflects 



thermal fluctuation in the critical region 
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. A peculiar feature of the A-shape of C(T) is 



that in addition to the sharp peak, the specific heat displays a symptom of its rise to the 



sharp peak already at T\ < T < T\ + 0.7K, which suggests the cause of the gradual fall of the 
shear viscosity in the same temperature region. A remarkable fact is that an ideal Bose gas 
with the same density as the liquid helium 4 shows a similar gradual rise of C(T) around T c 
as in Fig. 1(c). Normally, the enhanced thermal fluctuation at the vicinity of T c comes from 
strong competitions between the interaction energy and the entropy effect, and occurs only 
within the very narrow temperature region around T c . Hence, it is questionable to apply the 
concept of thermal fluctuation to this gradual rise of the specific heat at T\ < T < 2.8K, 
because of the following reasons: (1) it occurs in a far wider temperature region than the 
critical one, and (2) C{T) of an ideal Bose gas suggests that it really occurs even if there is 
no particle interactions. Since the gradual fall of i](T) above T\ may come from the same 
mechanism as the gradual rise of C(T) above T\, the thermal fluctuation is a questionable 
interpretation for both phenomena. (This interpretation of the A-shape of C(T) leads us to 
reconsider the nature of the A transition and the meaning of the infinite volume limit. See 
Sec.5A.) 

Compared to an ordinary liquid, a liquid helium 4 above T\ has a 10 -3 times smaller 
shear viscosity coefficient. Although in the normal phase, it is already an anomalous liquid 
under the strong influence of Bose statistics. Hence, it seems natural to assume that large 
but not yet macroscopic coherent wave functions gradually grows as a thermal equilibrium 
state as T — > T\ in the normal phase, and suppresses the shear viscosity [l^jisj]. (At an 
early stage in the study of superfluity, R.Bowers and K.Mendelssohn d.seussed a similar 
view on the nature of the decrease of r\ above T\ |3j.) The existence of such wave functions 
must affect the non-dissipative flow as well, such as the rotational flow above T\. Recently, 
from this viewpoint, the rotational property of a liquid helium 4 above T\ is theoretically 
examined [IjJ], and a slight decrease of the moment of inertia above T\ is predicted. (The 
similar phenomenon occurs in the Meissner effect of the charged Bose gas as well 15|.) 

Experimentally, although a lot of measurements have been done on superfluidity of a 
liquid helium 4, there is so little direct experimental information about flow patterns. Among 
various visualization techniques used in ordinary liquids, particle image velocimetry (PIV), 
which records the motion of micrometre-scale solid particles suspended in the fluid as tracer 
particles, recently becomes available in a superfluid helium 4 16j] 171 ] ■ Using this technique, 
we will be able to closely compare the theory and the experiment on the onset mechanism 
of superfluidity. 
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In this paper, after formulating the shear viscosity of a liquid, we generalize the result of 
Ref.14 to the capillary flow. Since we focus on the continuous change of the system around 
T\, we cannot assume from the beginning the sudden emergence of the macroscopic con- 
densate at T\. Rather, beginning with the system without the condensate, we consider the 
Bose system above and below T A on a common ground, and study the intricacy underlying 
the onset of superfluidity. 

This paper is organized as follows. Section 2. A considers the shear viscosity of a liquid 
by relating the dissipative flow to the non-dissipative ones using Kramers-kronig relation, 
and obtains a formula describing the influence of Bose statistics on the dissipative flow. 
Section. 2. B discusses the physical reason that Bose statistics suppresses the shear viscosity. 
Beginning with the Bose system without the condensate, Sec. 3 formulates a perturbation 
expansion of the susceptibility with respect to the repulsive interaction. By taking peculiar 
diagrams obeying Bose statistics in the above formula, Sec. 3 examines how the growth of 
the coherent wave function gradually suppresses the shear viscosity at T > T\. Section 
4 makes a comparison to experiments, and Sec. 5 discusses the nature of the A transition, 
compare the dissipative and the non-dissipative flows, and makes a brief comment on the 
shear viscosity of Fermi liquids. 



II. SHEAR VISCOSITY OF A BOSE LIQUID 
A. Formalism of the shear viscosity of a liquid 

For a liquid helium 4, the repulsive particle picture is not so unrealistic as it would be 
for any other liquid. Hence, as its simplest model, we use 

H = J2 <P)H% + K-^W^p'^ ( u > °)> ( 4 ) 

v p,p' <? 

where <3> p denotes an annihilation operator of a spinless boson. 

In the linear response theory, the shear viscosity coefficient is given by the following 
two-time correlation function of the stress tensor J xy (t) 

1 f°° 

V = J dt < J xy (t)J xy (0) > . (5) 

In a liquid, the particle interaction U enhances not only a few-particle or collective excita- 
tions but also their relaxations, hence reducing the relaxation time r, and the shear viscosity 



coefficient 77 = Gr. In principle, the shear viscosity coefficient of a liquid must be obtained 
by calculating an infinite series of the perturbation expansion of Eq. (5) [3] . The interaction 
U affects the imaginary part of the self energy of particles, thus changing the stress tensor 
J xy . To obtain the resulting change of < J X y{t)J X y{fy >, one must know the dissipation 
mechanism of a liquid. In contrast with the case of a gas, the mechanism responsible for 
the shear viscosity in a liquid is similar to the vacancy motion in a solid, and is therefore 
a highly inhomogeneous process in the microscopic scale. Hence, it is difficult to derive the 
shear viscosity of a liquid without using any phenomenological model 

Instead of 77, however, if we apply the linear-response theory to its reciprocal I/77, and 
make the perturbation expansion with respect to U, the increase of U naturally leads to the 
decrease of '77 , thus naturally describing the effect of U on 77. Along this line of thought, we 
regard Eq.(2) of the Poiseuille flow, in which 77 appears in the denominator, as a macroscopic 
linear-response formula. Without loss of generality, we may focus on a flow velocity at a 
single point. Using the flow velocity on the axis of rotational symmetry (z-axis), we define 
a current j = pv(r = 0) on the z-axis, and obtains 

where a is the conductivity of a liquid in a capillary flow [P = Pe z ). Equation. (6) is a 
longitudinal linear response of a liquid to AP/L, which includes the dissipation by the shear 
viscosity [^J. The perturbation expansion of o with respect to U properly incorporates to 
77 the strong-coupling nature of a liquid. 

Let us generalize Eq.(6) to the case of the oscillatory pressure gradient as follows 

3(u) = -o-(u)a 2 — (7) 

The Navier-Stokes equation gives an expression of the conductivity spectrum o~(uj) (see 
Appendix. B). By definition, a{uj) in Eq.(Y) must contain the mass m of particles. Specifically, 



o~{oj) satisfies the following sum rule [21] 



a{uj)duj = nm, (8) 







where n is the number density of particles. 

At T < T\, a frictionless flow appears in Eq.(7). Characteristic to the frictionless flow 
is the fact that in addition to the normal fluid part cr n (u), the conductivity spectrum cr(u) 

8 



p 

47? 



(a) 




CO (Wa 2 ) 



Act 



(b) 



A6(O0) 




CO 



FIG. 3: The change of the conductivity spectrum a(uj) from (a) at T > T\ to (b) at T < T\. At 
T > T\, a(u>) is given by Eq.(B5). At T < T\, u s is a half width of the sharp peak around uj = 
with an area A, and tu c is a frequency at the broad peak in a n (u). 

has a sharp peak around uj = (with a very small but finite half width u s ), while satisfying 
Eq.(8). Hence, one obtains 

,APM 



= - K(^) + A5(w)]a 2 



(9) 



where 5(u) is a simplified expression of the sharp peak around uj — 0, and A is an area 
of this peak [8|. Figure 3 schematically illustrates such a change of <j{oS) when the system 
passes T\. When a(u) has a form of (T n (c<J) + A5(co , ), one obtains the shear viscosity coefficient 
rj in Eq.(6) as 



V(T) = -- 



P 



A 



(10) 



<7 n (0) + 

Us 

In view of Eq.(10), the sharpness of the peak around u = (u s ~ 0) leads to the disap- 
pearance of the shear viscosity (r/(T) ~ 0) at T < T\. 

Let us generalize Eq.(7) so that it includes not only the current responding in phase 
with AF(w) but also a current responding out of phase. <j(uj) in Eq.(7) is generalized to a 
complex number ci + 102 as follows 

,AP(u) 



j(u) = - [<7i{uj) + i<7 2 {oj)\ a 



(11) 



where a{oj) in Eq.(7) is replaced by <J\{uS). Causality states in Eq.(ll) that only after the 
pressure gradient is applied, one observes a current, and obtains the following Kramers- 



Kronig relation 



a x {u') = - / duj 21 2 , 12 



a* u/ = / ^ 2 U ; /2 . 13 



We will embed Eq.(ll) in a more general picture, and derive properties of a capillary flow 
by regarding it as a special case. Instead of AP(w), we use a quantity directly specifying 
the state of the system, that is, the velocity field v(t). Applying the pressure gradient to a 
liquid is equivalent to assuming the velocity field v(t) satisfying the equation of motion 

dv(t) AP(t) 



dt 



(14) 



as an external field 
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Let us consider a rotational flow like Fig.2(b). We assume Vd(r,t) = ft x re wt of the 
rigid-body rotation. In the rotational flow, j(r,u) and Vd(r,u) form the perturbation 
energy like J dxVd ■ j. Since the rotational flow is a non-dissipative system, j(r,uj) is a 
dynamical response of a liquid to the mechanical external field Vd[r,Lj), such as j(r,u) = 
x(r,u)vd(r,u), where x( r *,^) is the susceptibility of the system [10(. Replacing AP(lu)/L 
with — iujv(uj) in Eq.(ll), we interchange the real and imaginary part of the coefficient in 
Eq.(ll), and make it a local equation as follows 

j{q,u) = \-u(y 2 (q,uj) + iuai{q,u)}a 2 v(q,uj). (15) 

Equation (15) is a general formula including not only iuai(q, uS) for the dissipative flow 
(Fig. 2(a)), but also — co , o 2 (g, uj) for the non-dissipative one (Fig. 2(b)). If one determines 
02(0, oj) in the non-dissipative flow, one obtains ai(q,u) using Eq.(12) at a given q. 

In the rotational flow, because of divv,i(r,t) = 0, v^r, t) acts as a transverse-vector 
probe to the excitation of bosons. This fact allows us a formal analogy with the vector 
potential A{r) in the Coulomb gauge (divA(r) = 0) acting on the charged Bose system. 
Hence, the neutral Bose system has a similar form of current to that of the charged Bose 

n 

system [1CJ such as 

p,n ^ 
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(h = 1 and t = it). Within the linear response, one obtains the generalized susceptibility 

1 fP 

XtwiQiVn) = — J dr exp(tuj n T) < G\T T J fl (q,T)J u (q,0)\G > , (17) 



where \G > is the ground sate 23|. Generally, the susceptibility is decomposed into the 



longitudinal and transverse part (/x = x, y, z) 

Xlxu (q,u) = ^fV(g,c) + (V - q -^\ X T (q,uj). (18) 
For the later use, we define a term proportional to q^q^ in XlMV by Xil/V 

'x L (g,uj) - x T (q,u) 

q 2 



= <W T (g,^) + Xi*v(q,v), (19) 



where XiW represents the balance between the longitudinal and transverse susceptibility 
241 ] . In the rotational flow, the influence of the wall motion propagates along the radial 
direction, which is perpendicular to the particle motion driven by rotation as illustrated in 
Fig.2(b). Hence, the rotational motion is a transverse response described by the transverse 
susceptibility x T (g, c<j). In the right-hand side of Eq.(15), the real part — a 2 co>o" 2 (<?, to) of the 
coefficient describes the rotational motion as a non-dissipative response such as 

- a 2 uja 2 (q,uj) = x T (q,u), (20) 

and the imaginary part a 2 uj<Ji(q, to) describes the dissipation. 

Let us focus on the capillary flow. Using Eq.(20) in the right-hand side of Eq.(12) at a 
given q, one obtains 

a = --J o dv-jr—z. (21) 

The slow excitations (u ~ 0) in x T {q, UJ ) of Eq.(21) are important for determining a\(q,u). 
On this point, Eq.(21) has different interpretations in the normal and the superfluid phase. 
In the normal fluid phase, x L (g,u;) = x T (g,co>) is satisfied at small q and u, and one can 
replace x T (g, tu) in Eq.(21) by x L {Qi 1jJ ) at a small q and to [25]. Hence, the conductivity 
(Jin(q, uj) of a capillary flow in the normal fluid phase is given by 

a a ln (q,u) = / duj— -. (22) 

7T /n UJ 1 — UJ' 1 
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In the superfluid phase, under the strong influence of Bose statistics, the condition of 
X L (q, of) = x T (q,uj) at q — > is violated in the low-energy region extending from uj =0 (see 
Sec.2.B). (Furthermore, the rigidity of superfluidity requires that the nonzero x L {q-i UJ ) ~ 
X T {q, u >) depends on q and uj very weakly at small q and a;.) Consequently, one cannot 
replace x T (?, <^) by x L {q-i UJ ) m Eq.(21). For ai(q,u), in addition to cr ln (q,uj), one must 
separately consider the contribution of x T (g, cj) — x L (<?)^)- For ai(0,u) of the capillary 
flow, one obtains 

( r, f n 2 [°° \im q ^ [x L (q,uj)-X T (q^)] m , 

(Ti(u) = (T ln (u) H ^ duj , (23) 

na z J uj 1 — uj' z 

where we replace x L (g, w) — x T (g, a;) by its limit at q — > 0. Since this limit depends on uj 

very weakly at small u, we take out it from the integral with the aid of 

f^ = ,M - < 24 > 

and obtains 

aiH = a ln (u) + lim[x L (g, 0) - X T (q, 0)]S(u>). (25) 

Comparing Eq.(25) with Eq.(9), one obtains an expected form of <j(uj) for the superfluid 
flow with A = 2 \im. q ^Q[x L {q, 0) — x T (<L 0)]/ (ira 2 ). If superfluidity was perfectly rigid to any 
perturbations, \im q ^ [x L (q, uj) — X T (<?> ^)] would have no frequency dependence. In reality, 
this limit weakly depends on uj in Eq.(23). Hence, although the peak is well approximated 
by S(u) in Eq.(25), it has a small but finite half width uj s as in Fig.3(b), which represents a 
degree of the rigidity of superfluidity. 

The change of <j(uj) from Fig.3(a) to 3(b) occurs under the sum rule Eq.(8). The area A of 
the sharp peak is equal to that of the shaded region in Fig.3(b). We approximate this shaded 
region by a triangle. (A broad peak of a n {uj) is located at uj c ). At uj = 0, the conductivity 
is enhanced by the sharp peak as A/uj s , but the original cr(0) changes to <t(0) — Act by the 
sum rule. We take into account the sum rule as —Aauj c /2 + A = 0, hence Act = 2A/uj c . 
After all, the conductivity at uj = becomes [cr(0) — 2A/uj c ] + A/uj s at T < T\. Using this 
form and the expression of A in Eq.(25), one obtains rj = p/(4a(0)) as 

V(T) = (?) z-n ^ • (26) 



4/ 2/1 2 ... . L ,. 



ain(0) + — hm[x L (g, 0) - X 1 (q, 0)] 

ix a L \uj s uj c ) <?^o 

Here, we define the mechanical superfluid density p s (T) (= \im q ^ [x L (q, 0) — x T (?, 0)] = 
lim^of*? 2 / \q^qu)]x^v) , which does not always agree with the conventional thermodynamical 
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superfluid density p s {T). (By "thermodynamical" , we imply the quantity that remains finite 
in the V — > oo limit.) Using u s <C oj c in Eq.(26), we obtain the following formula of the 

kinematic shear viscosity v{T) = r)(T)/p(T) 

v( T ) = /% m , (27) 

7ia z uj s p 

where v n is the kinematic shear viscosity in the normal fluid phase, and it satisfies <7i n (0) = 
l/(4i/ n ) in Eq.(26). 

One notes the following features in the formula (27). 

(1) v(T) of a superfluid is expressed as an infinite power series of u n , and the influence of 
Bose statistics appears in its coefficients. This result does not depend on the specific model 
of a liquid, but on the general argument (see Sec.5.B). (The microscopic derivation of u n is 
a subject of the liquid theory.) 

(2) Because of oo s ~ 0, which characterizes the sharp peak in Fig. 3(b), a small change 
of p s (T) in the denominator of the right-hand side is strongly enhanced to an observable 
change of v(T). 

(3) The existence of 1/a 2 in front of p s (T)/p indicates that a frictionless superfluid flow 
appears only in a narrow capillary with a small radius a: A narrower capillary shows a 
cleaner evidence of a frictionless flow. 

We can simply examine the condition of p s {T) ^ by ignoring the repulsive interaction. 
Using Eq.(16) in Eq.(17), one obtains the following term proportional to q^q v in x^v 



X "" M - ~V \ u + e (p)- e(p + (( ) • < 28 » 

where f{e{p)) is the Bose distribution. When bosons would form the condensate, f{e{p)) 
in Eq.(28) is a macroscopic number for p = and nearly zero for p ^ 0. Thus, in the sum 
over p in the right-hand side of Eq.(28), only two terms corresponding to p = and p = —q 
remain, with a result that 

X,,M)=mn c (T)^, (29) 

where mn c (T) is the thermodynamical superfluid density p s (T), and n c (T) is the number 
density of particles participating in the condensate. In view of Eq.(29), we find that the 
mechanical superfluid density p s (T) in this case is equal to the thermodynamical superfluid 
density p s (T). 

13 



When bosons form no condensate, the sum over p in Eq.(28) is carried out by replacing 
it with an integral, and one notes that q~ 2 dependence disappears in the result. Because 
of p s {T) = 0, one obtains v(T) = v n in Eq.(27). This means that, without the interaction 
between particles, BEC is the necessary condition for the superffuid flow. Equation. (27) 
shows that unless the particle interaction is taken into account, a superffuid flow abruptly 
appears at T\. To explain the gradual decrease of the shear viscosity just above T\ as in 
Fig. 1(a), we must obtain x T (g, cj) — x L (g, tu) under the repulsive interaction U. 



B. Bose statistics and repulsive interaction 

There is a physical reason to expect that the shear viscosity of a Bose liquid falls off at 
low temperature. Quantum mechanics states that in the decay from an excited state with 
an energy level £ to a ground state with E , the higher excitation energy E gives rise to 
the shorter relaxation time r. The time-dependent perturbation theory gives us 

-~\E-E \. (30) 
r 

In Eq.(21), <jx(q,u) = p/(4Gr) in the left-hand side includes the relaxation time r, whereas 
X T {q,uj) in the right-hand side includes the excitation spectrum, hence the difference of 
energy \E — E \. In this sense, Eq.(21) is a many-body theoretical expression of Eq.(30). 
Microscopically, the structural characteristic of liquids is the irregular arrangements of their 
molecules. These arrangements are similar, and therefore their energies differ only slightly 
from one arrangement to another. The transformation from one arrangement to another 
with a small \E — E \ continuously occurs with a large r. Since the fall of the shear viscosity 
at the vicinity of T\ is attributed to the decrease of r, Eq.(30) says that the excitation 
energy E must rise owing to Bose statistics. The relationship between the excitation energy 
and Bose statistics dates back to Feynman's argument on the scarcity of the low-energy 



excitation in a liquid helium 4 [26j], in which he explained how Bose statistics affects the 
many-body wave function in configuration space. To the shear viscosity, we will apply his 
explanation. 

Consider the velocity field in Fig. 2 (a) (depicted by long thin arrows) that moves white 
circles on a solid straight line to black circles on a one-point-dotted-line curve. (A viscid 
liquid would show such a spatial gradient of the velocity field. The influence of adjacent 
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layers in a flow propagates along the direction perpendicular to the particle motion. Hence, 
the excitation caused by the shear viscosity is a transverse excitation.) Let us assume that 
a liquid in Fig. 2 (a) is in the BEC phase, and the many-body wave function has permutation 
symmetry everywhere in a capillary. At first sight, these displacements by long arrows 
seem to be a large-scale configuration change, but this result is reproduced by a set of 
slight displacements by short thick arrows. In contrast with the longitudinal displacements, 
the transverse displacement does not change the particle density in the large scale, and 
therefore it is always possible to find, in the initial configuration, a particle being close to 
the particle after displacement [27|. In Bose statistics, owing to permutation symmetry, one 
cannot distinguish between two types of particles after displacement, one moved from the 
neighboring position by the short arrow, and the other moved from distant initial positions 
by the long arrow. Even if the displacement made by the long arrows is a large displacement 
in classical statistics, it is only a slight displacement by the short arrows in Bose statistics . 

To grasp the whole situation, it is useful to imagine this situation in the 3N-dimensional 
configuration space. The excited state made by short arrows lies in a small distance from the 
ground state in configuration space. Since the excited state is orthogonal to the ground state 
in the configuration integral, the many-body wave function corresponding to the excited state 
must spatially oscillate. Accordingly, it oscillates within a small distance in configuration 
space. Since the kinetic energy of the system is determined by the 3N-dimensional gradient 
of the many-body wave function, this steep rise and fall of the amplitude implies that the 
excitation energy is not small. 

In coordinate space, although the existence of irregular arrangements of molecules, whose 
energies only slightly differ from that of the ground state, is a structural characteristic of a 
liquid, the above result means that its number remarkably decreases due to Bose statistics, 
leaving only excited states with not low energies. The relaxation from such an excited 
state is a rapid stabilization process with a small r. This mechanism explains why Bose 
statistics reduces the shear viscosity coefficient rj = Gr 28| . Formally, while the longitudinal 
excitations are ensured by the particle conservation to satisfy x i (o , )0) = p in Eq.(25) both 
above and below T\, the low-energy transverse excitations become scarce below T\, and it 
reduces x T (<L 0) at q — > 0. Hence, o~i{u)) has a 5-function form in Eq.(25) and p s {T) becomes 
nonzero, thus leading to v{T) — > in Eq.(27). 

In Sec.2A, we considered the rotational flow as a complementary phenomenon of the 
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capillary flow. As discussed in the argument above Eq.(20), the particle motion by rotation 
is a transverse excitation. In Fig.2(b), the similar argument holds on the role of Bose 
statistics: Two types of particles, one from distant place by the long curved arrow, and 
the other from neighboring place by the short arrow are indistinguishable, hence leading to 
the scarcity of the low-energy transverse excitation as discussed in Ref.14. Under the slow 
rotation, the region near the center of rotation decouples from the motion of a container 
[3^]. A common mechanism lies in the two typical manifestations of superfluidity, the 
frictionless flow in a capillary and the decrease of the moment of inertia in the rotational 
flow in a container. 

When the system is at high temperature (/x <C 0), the coherent wave function has a 
microscopic size. If a long arrow in Fig. 2 takes a particle to a position beyond the coherent 
wave function including that particle, one cannot regard the particle after displacement as 
an equivalent of the initial one. The mechanism below T\ in Fig.2 does not work for the 
large displacement extending over two different wave functions. Hence, we obtain x T (<2S 0) = 
X L {q, 0) at q -> in Eq.(25), and u(T) = u n in Eq.(27). 

From our viewpoint, a remarkable state of the boson system lies at the vicinity of T\ 
in the normal phase (/i — S < 0). In Fig.4, the coherent many-body wave function grows 
to a large but not yet a macroscopic size. (Approximately, the chemical potential /i, hence 
fj, — E as well, is inversely proportional to the size of the coherent many-body wave function 
. which corresponds to the size of regions enclosed by a dotted line in Fig.4.) The 
permutation symmetry holds only within each of these regions. When particles are moved 
from a region A to another region A' in Fig.4(a), the mechanism below T\ does not work. 
But the repulsive interaction U, not only determining v n in Eq.(27), but also affects the role 
of Bose statistics as follows. In general, when a particle moves in the interacting system, 
it induces the motions of other particles. In particular, the long-distance displacement of a 
particle in coordinate space causes the excitation of many particles, and therefore it needs a 
high excitation energy. Hence, the short- distance displacement is a major ingredient of the 
low-energy excitation. This means that excited particles are not likely to go beyond a single 
coherent wave function, but likely to remain in it. Hence, under the repulsive interaction, 
the suppression mechanism of the low-energy excitation owing to Bose statistics works just 
above T\ as well. (A similar argument holds in the rotational flow in Fig.4(b) 14|.) 

This interplay between Bose statistics and the repulsive interaction will appear as follows. 
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FIG. 4: Two types of flows just above T\. Areas enclosed by a dotted line represent large but not 
yet macroscopic size of coherent wave functions, a size of which is exaggerated for clarity 

If we increase the strength of U, the excited bosons get to remain in the same coherent wave 
function, and therefore the low-energy transverse excitation will raises its energy. The 
relaxation time to the ground state becomes short, thus making the shear viscosity fall off. 
In Eq.(26), the condition of x L (<?> 0) = x T {Qi 0) at g ^ will be violated at a certain critical 
value of U . In real materials, we realize this situation in an alternative manner. If we 
decrease the temperature at a given U, the coherent wave function obeying Bose statistics 
grows. The above condition will be violated at a certain temperature , which will give the 
onset temperature of the decease of rj in Fig. 1(a). 



III. SUSCEPTIBILITY AND SHEAR VISCOSITY ABOVE T\ 

Let us formulate the influence of Bose statistics on the shear viscosity. Equation (26) 
says that x L (<?> ^) — X T (Qi UJ ) i s a crucial quantity for this purpose. The problem is how to 
formulate in this quantity the interplay between Bose statistics and the repulsive interaction. 
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FIG. 5: The process of line exchanges due to Bose statistics occurring in J^J u exp[— J (1tHj(t)] 
of Eq.(31). Two times of the line exchange in Fig.5(a), between a current-current response tensor 
JfiJu (a large bubble) and a bubble excitation by the repulsive interaction U (an inner small bubble 
with a dotted line), yields Fig. 5(b). (The black and white small circle represents a vector and scalar 
vertex, respectively.) 



Under the repulsive interaction Hi(t), one considers the integrand of Eq.(17) as 



<G\T t JJx,t)J v (0,0)\G>-- 



< 0\T T J tJi (x,T)J u (0,0)exp 


l-f 

. Jo 


dT Hj(t ) 


|0 > 


< 0\exp 


- / dTHj(T) 

. Jo 


|0 > 



(31) 



Figure 5(a) illustrates the current-current response tensor J^(x, t)J v (0, 0) (a large bubble 
with ji and v) in the medium: Owing to exp{— / Hj{r)dr) in Eq.(31), scatterings of particles 
frequently occur in \G > as illustrated by an inner small bubble with a dotted line U. (The 
black and white circle represents a coupling to v d (q) and U, respectively.) The state \G > 
includes many interaction bubbles like the small one with various momentums p' and p' + q'. 
A solid line with an arrow represents 

1 



G(iuj n ,p) = 



(32) 



iuj n - e(p) - E + pi 

(fi is a chemical potential implicitly determined by l /_1 S[exp( ( 5[e(g) + £ — /i]) — 1] _1 = n). 
Owing to the repulsive interaction, the boson has a self energy S (> 0) (we ignore its uo 
and p dependence by assuming it small). With decreasing temperature, the negative \i 
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at high temperature approaches a small positive value of E, finally reaching Bose-Einstein 
condensation satisfying p, — E. 

When the system is just above T A in the normal phase, particles in the ground state \G > 
experience the strong influence of Bose statistics. Hence, the perturbation of J^x, r)J I/ (0, 0) 
in Eq.(31) must be developed in such a way that, as the order of the perturbation in- 
creases, it gradually include a new effect due to Bose statistics. Specifically, the large 
bubble Jfj,{x, r) J„(0, 0) and the small bubble in Fig. 5(a) form a coherent wave function as 
a whole. When one of the two particles in the large bubble and in the small bubble have 
the same momentum (p = p'), and the other in both bubbles have another same momentum 
(p + q = p' + q') in Fig.5(a), a graph made by exchanging these particle lines must be in- 
cluded in the expansion. The exchange of two lines having p and p'{= p), and the exchange 
of the other two lines having p + q and p + q'{— p + q) yield Fig. 5(b). The result is that 
two bubbles with the same momentums are linked by the repulsive interaction, yielding the 
following contribution to Xfiv(q, w) 

*4z(p +§)/*(?+§)" 

V 

(a) With decreasing temperature, the coherent wave function grows to a large size, and 
the exchange of particle lines owing to Bose statistics like Fig. 5 occurs many times. Hence, 
one cannot ignore the higher-order terms in Eq.(31), which become more significant with 
the growth of the coherent wave function. Continuing these exchanges, one obtains the l-th. 
order term which has two vector vertex (black circles) at both ends, to which (p+ §)^G»+ f)^ 
is attached, and 21 scalar vertex (white circles) between the two ends. 

(b) Among physical processes in Eq.(33), it becomes evident with decreasing temperature 
that a process including zero momentum particles plays a dominant role. Specifically, a 
process with p = corresponds to an excitation from the rest particle, and a process with 
p = —q corresponds to a decay into the rest one. These processes are important in the 
higher-order terms as well. 

With these considerations (a) and (b), we obtain the following form of Xfiv{q,Q) = 
[x L (q, 0) - x T (q, ^)}q^qu/q 2 at the vicinity of T A 

1 oo 

v(?.o) = ^E^(«) w ' ( 34 ) 

1=0 
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/( £ (p) + S)-/(e(p + g ) + £ 
oo + t{p) - tip + q) 



z 



(33) 



where 



(exp(/3[S - /i] - I))' 1 - (exp(/?[e(g) + g -//])- I) 



(35) 



is a positive monotonously decreasing function of q 2 , which approaches zero as q 2 — > oo. 

At a high temperature <C 0), Fp(q) is small, and it guarantees the convergence of an 
infinite series in Eq.(34), with a result that 

QiiQu 1 Fp(q) 



(36) 



2 Vl-UF p { q y 

With decreasing temperature, however, the negative p gradually approaches S, hence S — 
/i — > 0. Since Fp(q) increases as S — /i — > 0, it makes the higher-order term significant in 
Eq.(34). Since is a positive decreasing function of q 2 , the divergence of Eq.(36) first 

occurs at q 2 = 0. An expansion of Fp(q) around q 2 = 0, i*]a(0) — 6g 2 H — ■ has a form such as 



Fp(q) 



P 



4smh2 /WT)-E] 



1 - 







1 



2 ta,h('^ (r 2 ) - S ^ 2 '" 



+ 



(37) 



At q — > 0, the denominator 1 — UFp(q) in the right-hand side of Eq.(36) has a form of 
[1 — UFp(0)] + Ubq 2 . With decreasing temperature (S — — > 0), UFp(0) increases and 
finally reaches 1, that is, 



UP = 4 sinh 



2 //5[/i(T) - E(U)} 



(38) 



V 2 

At this point, the denominator in Eq.(36) gets to begin with q 2 , and Xuu(q,0) there- 
fore changes to have a form of q^q v /q 2 at q — > 0, and q^q v /q 2 has a non-zero coefficient 
F (3 (0)/(2V r f/6). In view of the definition of Eq.(19), this result implies lim q ^ [x L (q, 0) — 
y T (q, 0)] in Eq.(26). From now, we call T satisfying Eq.(38) the onset temperature T on 



|3l|. 



At T = T on , substituting Eq.(37) into Eq.(36), we obtain x^v at q — > 



2,/ ' 1 tanh ^ 1^°" ^^°") ~ 9/*?" 



r t,r V 2 y r 

By the definition of p s (T) = \im q _, (q 2 / q^q^x^v, one obtains with the aid of Eq.(38) 

1 



(39) 



rn 



^sinh \P on [p(T on ) - E] 



(40) 
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The mechanical superfluid density p s (T) is given by mc(T)no(T), where 

exp([3\p(T) -E|) + 1 
is a Fermi-distribution-like coefficient, and 

"» (T) = Fexp(- ;3W r)-E])-i (42) 

is the number density of p = bosons. In the V — > oo limit, this quantity is normally 
regarded to be zero. For the finite system just above T\, however, no(T) has a large but not 
yet macroscopic value. In real finite system, its magnitude must be estimated by experiments 
(see Sec. 4). 

The violation of x T (<7> 0) = X L {Qi 0) at g ^ gives rise to various anomalous mechanical 
responses of a superfluid. At the vicinity of T\, Eq.(38) is approximated as U(5 = (3 2 [p{T) — 
E([/)] 2 for a small p — E. This condition has two solutions p{T) = E(C7) ± y/TTk^T. It 
is generally assumed that the repulsive Bose system undergoes BEC as well as a free Bose 
gas. Hence, with decreasing temperature, p(T) of repulsive Bose system should reach E(£7) 
at a finite temperature, during which course the system necessarily passes a state satisfying 
p(T) = E(C/) — y/UksT '. Consequently, the anomalous mechanical response of a superfluid 
appears prior to T\. Using Eq.(40) in the formula of the kinematic shear viscosity (Eq.(27)), 
we obtain v(T) at the vicinity of T\, which always falls off prior to BEC , that is, T on > T\. 
Among various mechanical responses of a liquid helium 4 above T\, the shear viscosity 
coefficient most evidently shows a change toward superfluidity far above the critical region. 
This feature comes from the structure of Eq.(27), in which the small us s amplifies a small 
p s (T)/p to an observable change of f(T). 

The phenomenon at T on has the following difference from thermal fluctuations. Since the 
thermal fluctuation is essentially a local phenomenon, its influence already appears in the 
low-order terms of the perturbation expansion in which only a few particles participate. In 
Eq.(34), however, only after the perturbation expansion is summed up to the infinite order, 



the non-zero p s (T) defined as \im g ^ [q 2 / (q^q u )]x^ u appears [12| [32 ]. 

With decreasing temperature from T on , in addition to the particle with p = 0, other 
particles having small but finite momentums get to contribute to the 1/q 2 divergence of 
Xhu{q,0) as we ll- (I n addition to Eq.(35), a new Fp(q) including p ^ also satisfies 1 — 
UFp{Qi) = in Eq.(36).) Hence, the mechanical superfluid density increases. On this change, 
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we have the following physical explanation. For the repulsive Bose system, particles are likely 
to spread uniformly in coordinate space due to the repulsive force. This feature makes the 
particles with p ^ behave similarly with other particles, especially with the particle having 
zero momentum. If they behave differently from others, a resulting locally high density of 
particle raises the interaction energy. This is a reason why many particles participate in the 
singular dynamical behavior even at a temperature in which only a few particles participate 
in the coherent many-body wave function . In contrast with the thermodynamical superffuid 
density p s {T), the mechanical superffuid density p s (T) is a concept including such a dynamics 
of the system. 

When the system reaches T = T\ in which p, = £ is satisfied, one notes c(T\) = 1 
and 7io(T\) = n c . This means that the zero momentum part of p s {T\) agrees with the 
conventional thermodynamical superfluid density p s {T\) = mn c , and n c abruptly grows to 
a macroscopic number. While p s {T) at T > T\ vanishes in the V — > oo limit, p s (T) at 
T < T\ remains finite in this limit. Hence, in thermodynamics, only the latter remains. 
In the mechanical response of the finite system, however, p s {T) manifests itself already at 
T > T\. 



IV. COMPARISON TO EXPERIMENTS 



At T > 2.8K in 1 atm, the shear viscosity of a liquid helium 4 slightly increases with 
decreasing temperature. As shown in Fig.l(a), after reaching a maximum value at 2.8K, 
it begins to reduce its value. On the contrary, the shear viscosity of a classical liquid 
has the general property of increasing monotonously with decreasing temperature (33]. At 
first sight, a liquid helium 4 seems to undergo a crossover from a classical to a quantum 
liquid at 2.8K, but it is superficial. Comparing Eq.(26) and (27), one notes that v(T) 
is a more appropriate quantity than r](T) for describing simply the change of the system 
around T\. As shown in Fig. 6 (a), the total density p(T) of a liquid helium 4 in the normal 



phase monotonously increases with decreasing temperature until T = T\ [34|. Fig. 6(b) 
shows the kinematic shear viscosity v(T) = n{T) / p{T) using n{T) (Fig.l(a)) and p{T) 
(Fig.6(a)), revealing that v(T) monotonously decreases with decreasing temperature over 
the whole range of T\ < T < A.2K. This means that in 1 atm, it is just below the gas-liquid 
condensation point T c = 4.2i^ that the strong influence of Bose statistics begins to suppress 
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FIG. 6: (a) The temperature dependence of the total density p of a liquid helium 4 (after the data 
in Ref.34). (b) The temperature dependence of the kinematic shear viscosity v of a liquid helium 4 
obtained using Fig.l(a) and Fig. 6(a). (c) p s (T)/p obtained by Eq.(43) using v{T) in Fig. 6(b) and 
uj s = 5 x 10~ 3 rad/s. 

the shear viscosity in a liquid helium 4. 

This temperature dependence and the smallness of v{T) implies that a liquid helium 4 in 
the normal phase is already an anomalous liquid under the strong influence of Bose statistics. 
In a strict sense, it is not certain whether there is a completely classical regime in liquid 
helium 4 at 1 atm. But, even if we assume T on > 4.2K, it is impossible to analyze the onset 
mechanism of the decrease of v in a virtual liquid at T > A.2K in 1 atm. For the present, 
we assume T on = 4.2K, and interpret v(T) in Fig. 6(b) using Eq.(27) with v n = v{T on ). One 
obtains p s (T)/p by 

u{T)- l -u{T on )-\ (43) 



Ps(T) 



ira^uOg p 

A typical value of the capillary radius used in the experiment in Ref.3, 4 and 5 is a ~ 
5 x 10 _3 cm. In the rotation experiment by Hess and Fairbank 29|, the moment of inertia 
I z just above T\ is slightly smaller than the normal phase value If (see Sec.5.B). Using 
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these data, Ref.14 roughly estimates p s (T x + 0.03K)/p = 8 x KT 5 , and p„(T x + 0.28K)/p = 
3 x 10~ 5 . With these values and v(T) in Fig.6(b), we obtain a rough estimation of u s as 
w s ~5x 10~ 3 rad/s at these two temperatures just above T\. Figure 6(c) shows a resulting 
temperature dependence of p s (T)/p. (Since the precision of these two values of p s (T)/p 
derived from currently available data is limited, the absolute value of p s (T)/p in Fig.6(c) 
has a statistical uncertainty.) 

The quantity which directly indicates the onset of superfluidity in dissipative system 
is the change of the conductivity spectrum a(u). In view of the gradual decrease of 77 
above T\, the sharp peak around uj = and the corresponding change from a{uj) to a n (ou) 
in Fig. 3 must already appear at T > T\. (Since this change does not occur in the case of 
thermal fluctuations, it is useful for ruling out the possibility of fluctuations.) To confirm this 
prediction, the time-resolved measurement of the oscillating flow velocity v(r,t) is necessary 
under the slowly oscillating pressure gradient. Such an experiment must be performed in a 
thin capillary with an inner radius of ~ 10 _1 mm. Recently, the PIV technique recorded 
the velocity of tracer particles in a flow of liquid helium 4, and mapped out the velocity field 
v(r) in a pipe [3] (it| . If the PIV experiment is performed under the above conditions, 
the change of cr(u) in Fig.3 will be observed. There is now ambiguity in interpretation of 
the movement of tracer particles If tracer particles interact only with the normal fluid 
flow and trace its velocity, the change from <j(oj) to cr n (co>) as in Fig.3 will be observed. The 
quantity A in Eq.(lO) is equal to the area of the shaded region a{ui) —a n (uj) in Fig.3(b). On 
the other hand, if tracer particles interact with the superfluid flow as well, the emergence 
of the sharp peak around to = as in Fig.(b) will be detected. In any case, using thus 
obtained A, one can determine u s by interpreting the experimental data of t](T) with Eq.(lO). 
Furthermore, using A = 2p s (T)/(ira 2 ), one will obtain an independent estimate of p s (T)/p 
at T > T\. Although such an experiment may be a difficult one, it is an experiment worth 
attempting. 
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V. DISCUSSION 



A. Interpretation of the A-shape of specific heat 

Bose-Einstein condensation differs from ordinary phase transitions in that it occurs with- 
out interactions between particles. At the early stages in the history of the study of phase- 
transition in the early part of twentieth century, almost all phase transitions were thought to 
have a A-like-shaped temperature dependence of the specific heat C(T) at the vicinity of T c . 
As the precision of experiments was improved, however, it became clear that C(T) in most 
phase transitions has, not the A-like shape, but a 5-function-like one. In a liquid helium 4, 
however, the precise measurements revealed that C(T) intrinsically has the A-like temper- 
ature dependence. The principal peak at T\ has been a subject of intensive studies of the 
critical phenomena. Furthermore, the gradual rise of C(T) above T A has been interpreted 
as a sign of enhanced thermal fluctuations specific to the Bose systems near T A . 

There is a reason that such a fluctuation-based interpretation was accepted on the gradual 
rise of C(T). In theories of the Bose system, the V — > oo limit is normally used. In this 
limit, the number density of particles with zero momentum is exactly zero at T > T\, and is 
finite only at T < T A . (In the equation of states N/V = g 3 /2(e^) + V^e^^l - e^) where 
g a (x) = J2 n x n /n a , the zero momentum part V~ x e^j{\ — e^ M ) has a finite value at V — > oo 
limit only when e /3/1 /(l — e^ 11 ) diverges, that is, \x = 0.) The V — > oo limit is necessary 
for the definition of the phase transition that the thermal average of the order parameter is 
exactly zero above T c . 

The V — > oo limit is a good approximation only when one can clearly distinguish between 
microscopic and macroscopic phenomena. In the Bose system, however, the distinction be- 
tween "microscopic" and "macroscopic" is not so obvious as in other phenomena. Since 
particles with zero momentum plays an important role at low temperature, the coherent 
many-body wave function with a macroscopic wave length plays a significant role in micro- 
scopic descriptions of the system. In this sense, "microscopic" and "macroscopic" coexist 
from the beginning. But, if we naively apply this V — > oo limit to the anomalous response 
at T > T\, the thermal average of the output quantity is set to be zero by V — > oo, and we 
are forced to regard all of them as being caused by thermal fluctuations. We must carefully 
examine the validity of regarding a whole sample as occupying an infinite volume. In the 
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thermodynamical viewpoint using V — > oo, the intermediate-sized coherent wave function 
plays a minor role above T\, but for the Bose system above T\, it substantially affects the 
system. Hence, we must experimentally estimate the magnitude of effects hidden by this 
7 — > oo limit. In this paper, we viewed the gradual decrease of rj above T\ from this point 



36j , and showed that the growth of the intermediate- sized coherent wave function gives rise 



to the decrease ofrj . 



B. Superfluidity in the dissipative and the non-dissipative systems 

As an example of the non-dissipative systems, we considered a rotational flow. In the 
rotational flow, the quantity directly indicating the onset of superfluidity is the moment of 
inertia 

- ?f) , (44) 

where if is its classical value 14j. On the onset mechanism of superfluidity, one can see 
the physical difference between the dissipative and the non-dissipative systems in Eq.(27) 
and (44). 

(a) In Eq.(44), p s (T) appears as a correction to the coefficient of the linear term of If. 
On the other hand, in the expansion of Eq.(27) with respect to u n , p s (T) appears in the non- 
linear higher-order terms of v n . Furthermore, with decreasing temperature, the higher-order 
terms become dominant in Eq.(27). 

(b) In Eq.(44), the change of the transverse excitations appears in I z without being 
enhanced, and therefore the nonzero p s {T) only slightly affects I z above T\. On the other 
hand, in Eq.(27), the change of the transverse excitations does not directly appear in u, 
but through the dispersion integral as in Eq.(23), which is a characteristic feature of the 
dissipative systems. During this process, a sign of superfluidity in the dissipative systems is 
enhanced to an observable scale. 

(c) In Eq.(27), the peak width u s reflects a degree of the rigidity of superfluidity, which 
is derived from the to dependence of ]im. q ^. \x L (q, u ) ~ X T (<?; w )]) whereas ui s does not appear 
in Eq.(44). This feature is characteristic to superfluidity in the dissipative systems, which 
manifests itself while showing the rigidity against the dissipation. The microscopic derivation 
of u s is a future problem. 

In principle, it is possible to begin with Eq.(5), and obtain Eq.(27) after perturbation 
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calculations. But it may require a knowledge of the dissipation mechanism in a liquid. In 
other words, the correct model of the dissipation mechanism of a liquid must reproduce 
Eq.(27) on the influence of superfluidity, because Eq.(27) is based on the general argument. 
In this sense, one can use this requirement as a criterion of the correct model of a liquid. 

As another example of superfluidity in the dissipative system, one knows the anomalous 
thermal conductivity of a liquid helium 4 at T < T\. Heat flow q is expressed as q = — kVT, 
where k is the thermal conductivity coefficient. At T = T\, k jumps to a 10 7 times larger 
value than k just above T\. At T\ < T < 2.8K, however, the gradual rise of k corresponding 
to the gradual fall of rj is not observed {3^]. Formally, we do not know a non-dissipative 
phenomenon complementary to the heat conduction, like the rotational flow in the shear 
viscosity. Hence, we can not apply the formalism in this paper to the onset mechanism of 
the anomalous thermal conductivity. This difference between heat conduction and shear 
viscosity on a formal level may have some implication on the absence of the gradual rise of 
k at T x < T < 2.8K. 



C. Comparison to Fermi liquids 

In a liquid helium 3, the fall of the shear viscosity at T c is known as a parallel phenomenon 
to that of a liquid helium 4. The formalism in Sec. 2 is also applicable to the shear viscosity of 
a liquid helium 3. For the behavior above T c , however, there is a striking difference between 
a liquid helium 3 and 4. The phenomenon occurring in fermions at the vicinity of T c is not a 
gradual growth of the coherent wave function, but a formation of the Cooper pairs from two 
fermions. (This difference evidently appears in the temperature dependence of the specific 
heat: C(T) of a liquid helium 3 just above T c does not show the symptom of its rise to the 
sharp peak.) Once the Cooper pairs are formed, they are composite bosons situated at low 
temperature and high density, and immediately jumps to the ESP or BW state. Hence, the 
shear viscosity of a liquid helium 3 shows an abrupt drop at T c without a gradual fall above 
T 

. n electron superconductivity, the fluctuation-enhanced conductivity a' is observed above 



T, 



0. 



In bulk superconductors, however, the ratio of a 1 to the normal conductivity o n is 
about 10 -5 at the critical region, and zero outside of this region. Practically, it is unlikely 
that thermal fluctuations create a large change of a at temperatures outside of the critical 
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region. 



APPENDIX A: MAXWELL'S RELATION 

Consider the shear transformation of a solid and of a liquid. In a solid, the shear stress 
F xy is proportional to a shear angle as F xy = G<f), where G is the modulus of rigidity. 
The value of G is determined by dynamical processes in which vacancies in a solid move to 
neighboring positions over the energy barriers. As increases, F xy increases as follows, 

dt dt y ' 

In a liquid, the flow motion rearranges the relative position of particles, reducing the shear 
stress F xy to a certain value. Presumably, the rate of such a relaxation is proportional to 
the magnitude of F xy , and one obtains 

dFgyy _ _ Fxy ( A2) 

dt dt t { ' 

In the stationary flow after relaxation, F xy remains constant, and one obtains 

G— = (A3) 
dt r v 1 

Figure. 7 shows two particles 1 and 2, each of which starts at (x, y) and (x, y+Ay) simultane- 
ously and moving along the x-direction. Consider a velocity gradient v x (y) along y direction. 
After At has passed, they (1' and 2') are at a distance of Av x At along the re-direction. As a 
result, the shear angle increases from zero to A<f>, which satisfies Av x At = AyAtfi as depicted 
in Fig. 7. Hence, we obtain 

dy dt y ! 

Using Eq.(A4) in Eq.(A3), and comparing it with Eq.(l), one obtains rj = Gt (Maxwell's 

relation). 

APPENDIX B: CONDUCTIVITY SPECTRUM a(u) 

The Navier-Stokes equation under the oscillating pressure gradient APe^/L is written 

in the cylindrical polar coordinate as follows 

dv _ ( d d \ APe iwt 
dt \dr 2 rdr J pL 

28 




FIG. 7: In a liquid flowing along the x-direction, owing to the velocity gradient along the y- 
direction, a small part of a liquid with a rectangular shape transforms to a shape whose cross 
section is a parallelogram. 



The velocity field has the following form 

APe iuJt 

with the boundary condition of v(a,t) = 0. Av(r,t) satisfies 



(B2) 



dAv(r, t) 
dt 



= v 



d d , s 



(B3) 



and therefore Av(r,t) has a form of Bessel function J (i\r) with A = (1 + i)y/u/(2v). 
Hence, 

APe iu * / J (iXr)\ , . 



iwpL \_ J (i\a) J 
At r = 0, the conductivity cr(u;) satisfying pv(0,t) = a(uj)a 2 APe lujt / L is given by 

/ \ 



1 



1 - 



(B5) 



Re cr(u;) is schematically illustrated in Fig.3(a). (Re cr(0) in Eq.(B5) agrees with Eq.(6). 
Im cr(u) gives an expression of cr 2 (uj) in Eq.(ll), but it does not agree with 02(0,0;) in 
Eq.(20), because it is derived from the phenomenological equation with dissipation like the 
Navier-Stokes equation.) 
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